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FOREWORD 

The  results  of  a  study  conducted  under  DA  Project  1P014501B14A05 , 

AMCMS  Code  5011.11.85300.04,  are  presented  in  thi°  report. 

The  design  of  a  weapon  system  provides  a  natural  setting  for  an 
optimization  problem.  The  design  requirements  stipulate  that  the  system 
is  to  perform  some  task  at  some  index  of  performance.  The  optimizer  then 
is  to  search  for  the  design  parameters  such  that  the  weapon  system  not 
only  performs  its  task,  but  also  maximizes  its  performance.  The  objective 
of  this  study  is  to  apply  a  relatively  new  s teepee t-descent  numerical  pro¬ 
cedure  to  an  artillery  design  problem  which  involves  the  dynamic  behavior 
of  a  105mm  howitzer  which  is  fired  while  resting  on  rubber  tires.  The 
tires  act  like  a  spring  during  the  firing  cycl-  which  causes  the  weapon 
to  leave  the  ground  so  that  the  likelihood  of  it  being  zeroed  in  for  the 
next  round  has  been  reduced  considerably.  The  purpose,  then,  will  be  to 
minimize  the  pitch  motion  of  the  weapon  by  obtaining  a  set  of  design  para¬ 
meters  which  are  subject  to  equality  as  well  as  inequality  constraints 
prescribed  by  design  requirements. 
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SUMMARY 

An  artillery  weapon  mounted  on  tire;  or  tracks  has  some  undesirable 
features  at  high  angle  fire.  Unlike  the  hard  mount  (weapon  rests  on  a 
base  plate)  the  flexible  mount  will  have  a  pitch  motion.  That  is,  during 
the  recoil  stroke,  the  tires  load  up  or  compress  and  act  like  a  spring, 
and  just  before  counter recoil  begins,  the  tires  begin  to  unload  sending 
the  weapon  off  the  ground.  Such  a  phenomenon  is  known  as  a  secondary 
recoil  effect.  The  control  rod  design  becomes  much  more  difficult  with 
this  secondary  recaii  effect  because  an  additional  acceleration  term 
enters  into  the  recoil  equations.  Also,  it  is  obvious  that  when  the 
weapon  comes  to  rest  the  likelihood  of  it  being  zeroed  m  for  the  next 
round  has  been  reduced  considerably.  The  purpose  of  this  study  is  to 
reduce  the  pitch  motion  of  the  weapon  and  at  the  same  time  determine  the 
orifice  areas  tor  the  control  rod  design. 

This  study  was  performed  on  the  XM164,  a  light  weight,  105mm  howitzer. 
The  presen.,  control  rod  design  for  short  recoii  (75  degrees  elevation! 
fields  approximately  six  inches  of  ,:hop".  Results  from  this  study  show 
that  between  45  and  86  per  cent  reduction  in  the  pitch  motion  is  possible 
(depending  upon  which  design  option  is  used)  by  determining  the  optimal 
shape  rod  force.  Once  this  rod  force  has  been  round,  the  orifice  areas 
can  be  determined. 

A  steepest-descent  numerical  procedure  will  be  used  to  minimize  the 
pitch  motion  of  the  weapon  along  with  satisfying  certain  design  constraints 
imposed  upon  the  system.  This  technique  starts  with  an  estimated  design, 


analyzes  it,  and  then  improves  on  the  design.  It  is  an  iterative  process 
and  at  each  iteration  an  improvement  is  made  until  no  significant  gains 
can  be  achieved. 

The  results  of  this  study  clearly  indicate  that  weapon  performance 
can  be  improved  by  using  methods  of  optimal  design. 


I.  Introduction 

Weapon  systems  of  today  and  of  the  future  are  becoming  more  complex 
and,  as  a  result  of  this  complexity,  the  engineer's  intuition  and  experi¬ 
ence  become  increasingly  more  difficult  to  apply  because  of  the  possible 
trade-offs  in  the  design  parameters.  Because  the  task  of  the  engineer 
becomes  more  difficult  in  meeting  requested  design  requirements  as  weapon 
systems  become  more  complex,  it  is  important  that  the  design  procedure 
be  represented  by  mathematical  modeling,  i.e.,  a  translation  of  the  physi¬ 
cal  description  of  the  problem  into  mathematical  terms.  Although  a 
mathematical  model  may  be  formulated,  the  solution  may  still  be  difficult 
to  obtain  for  several  reasons  The  model  itself  may  become  very  complex 
and  that  which  is  even  more  difficult  to  cope  with  is  the  fact  that  some 
of  the  parameters  may  only  be  engineering  estimates  based  on  past  experi¬ 
ence  or  perhaps  very  little  is  known  about  the  dynamic  behavior  of  a 
parameter.  Also,  the  solution  must  be  a  physical  realization  of  the 
mathematical  design.  In  short,  the  conversion  of  mathematical  theory 
into  an  engineering  accomplishment  may  not  be  an  easy  task. 

The  design  of  a  weapon  system  provides  a  natural  setting  for  an 
optimization  problem  assuming  a  knowledge  of  all  environmental  factors 
which  influence  the  design  process.  The  design  requirements  specify 
that  the  system  is  to  perform  some  task  at  some  index  of  performance. 

To  determine  the  optimum  solution,  the  concept  of  index  of  performance 
is  introduced  and  will  be  defined  as  the  functional  relationship  among 
the  system  characteristics.  The  optimizer  then  is  to  search  for  the 
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admissible  parameters  such  chat  tne  weapon  system  not  only  performs  its 
task,  but  also  maximizes  its  performance.  As  design  specifications  tend 
to  tighten,  it  becomes  increasingly  important  to  design  optimum  systems 
relative  to  some  performance  criterion,  and,  in  fact,  specify  that  per¬ 
formance  be  optimized. 

It  is  only  natural  then  that  the  methods  used  in  the  design  of 
optimum  systems  be  of  interest  for  these  are  the  analytical  tools  which 
will  determine  the  results  for  the  optimal  design  problem.  Because  of 
the  computer,  many  different  disciplines  have  provided  revolutionary  aids 
with  respect  co  analytical  cools  for  the  solucion  to  problems  that  were 
seemingly  hopeless  only  several  years  ago.  The  objective  of  this  study 
is  to  apply  the  relatively  new  technology  to  an  artillery  design  problem 
and  to  develop  a  method  which  will  aid  the  engineer  in  obcalning  design 
parameters  subject  to  certain  constraints  and  require  chat  the  perform¬ 
ance  of  the  weapon  be  optimal  in  some  sense. 


II.  Statement  of  the  Problem 

An  artillery  weapon  mounted  on  tires  or  tracks  has  some  undesirable 
features.  Unlike  the  hard  mount  (weapon  rests  on  a  base  plate),  the 
flexible  mount  will  have  a  pitch  motion..  During  the  recoil  stroke,  when 
the  weapon  is  fired  at  75  degrees  elevation,  the  tires  load  up  or  com¬ 
press;  and  when  counter recoil  begins,  the  tires  act  like  a  spring  and 
unload  sending  the  tirt.3  off  the  ground.  It  is  quite  obvious  that,  when 
the  weapon  comes  to  rest,  the  likelihood  of  it  being  zeroed  in  for  the 
next  round  has  been  reduced  considerably,  especially  for  high  rate  of 
fire  weapons.  This  phenomenon  is  known  as  a  secondary  recoil  effect 
because  an  additional  acceleration  term  enters  into  the  recoil  equations. 
Because  of  this  secondary  recoil  effect,  the  control  rod  design  becomes 
much  more  difficult.  For  short  recoil,  the  orifice  areas  in  the  control 
rod  are  designed  at  maximum  elevation  (75  degrees);  therefore,  when  ele¬ 
vation  is  mentioned  throughout  the  remainder  of  this  report,  it  refers 
to  maximum  elevation.  The  weapon  positioned  for  high-angle  fire  is 
shown  in  Figure  1. 

The  purpose  of  this  study  will  be  to  develop  a  systematic  control 
rod  design  procedure  characterized  by  mathematical  modeling  for  the  high¬ 
speed  digital  computer.  Conceptually,  it  will  be  one  phase  of  a  study 
that  will  give  the  optimal  weapon  which  meets  the  given  design  require¬ 
ments.  To  do  this,  a  steepest-descent  numerical  procedure  will  be  used 
to  minimize  the  hop  or  pitch  motion  of  the  weapon  and,  at  the  same  time, 
to  determine  the  necessary  control  rod  design  which  wiii  minimize  hop. 
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A  second  phase  will  be  to  incorporate  geometrical  effects  into  the  optimal 
design  problem  in  order  to  establish  optimal  geometries  for  certain  con¬ 
figurations. 

The  recoil  equation  is  of  the  form  of  Equation  (1) 

x  +  f(x)x2  -r  g(x)  =  h(t)  (1) 

for  a  rigid  irmnt.  In  the  second  term  of  this  equation,  the  expression 
for  the  effect  of  the  control  rod  orifice  area  is  defined;  however,  with¬ 
out  any  loss  of  generality,  the  control  rod  orifice  areas  can  also  be 
obtained  from  a  predetermined  rod-pull  force  R(t).  For  the  flexible 
mount,  the  above  equation  is  coupled  with  the  equation  describing  the 
pitch  motion  of  the  weapon  and  thus  yielding  two  second-order  nonlinear 
ordinary  differential  equations  with  prescribed  initial  conditions.  The 
orifice  areas  are  a  function  of  the  state  of  the  system.  To  eliminate 
state  variable  inequality  constraints , .R(t)  will  be  taken  as  the  control 
variable  which  is  to  be  determined  to  minimize  hop  (the  pitch  motion  of 
the  weapon)  subject  to  other  design  constraints. 

This  study  was  performed  on  an  existing  weapon,  namely,  the  XM164. 

The  XM164  is  a  lightweight,  split-trailed  towed  105mm  howitzer  with  the 
XM44  hydropneumatic  recoil  mechanism.  Unlike  a  rigid  mount,  the  XM164 
is  a  flexible  mount  and  is  fired  while  resting  on  rubber  tires.  For  a 
rigid  mount  weapon,  the  resisting  force  R(t)  on  the  recoiling  parts  is 
designed  with  a  trapezoidal  shape  as  shown  in  Figure  2.  With  the  proper 
design  of  the  control  rod  orifice  area,  the  flow  of  oil  in  the  recoil 
mechanism  is  controlled  and  such  a  force,  as  shown  m  Figure  2,  can  bp 
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Figure  2.  Rod  Forte  During  Recoil  for  a  Rigid  Mount 
obtained.  However,  when  a  force  tshapea  as  m  Figure  2)  is  designed  for 
the  flexible  mount,  the  question  is  asxed,  "Can  this  force  be  applied 
with  some  other  'best-'  shape  such  that  it  will  reduce  the  pitch  of  the 
weapon?"  This  is  the  basic  question  with  which  this  study  is  concerned. 

In  this  report,  the  optimum  rod  force  is  defined  as  that  curve 
which,  according  to  some  measure  vthe  ncp  motion),  satisfies  all  of  the 
requirements  imposed  upon  the  system 
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III.  Formulation  of  the  Problem 

During  the  recoil,  counter recoil  cycle  there  are  four  different 
tines  which  are  of  concern.  These  are  shown  in  Figure  3.  At  these 
four  times 


Figure  3 

* 

t  -  initial  conditions 
0 

tj  -  end  of  the  recoil  stroke 

t  -  time  at  which  maximum  hop  occurs 
2 

t,  -  end  of  counterrecoil 
f 

certain  conditions  must  be  satisfied  from  the  design  requirements.  At 
time  tQ  the  initial  conditions  for  the  state  of  the  system  are  given. 

At  time  tj  the  displacement  of  the  recoiling  parts  is  required  to  be 
equal  to  some  specified  value  and  the  velocity  of  the  recoiling  parts 
must  be  equal  to  zero.  At  time  t^  the  velocity  of  the  pitch  motion  must 
be  zero  and  the  displacement  of  the  pitch  motion  is  to  be  a  minimum. 

Note  that  it  will  be  possible  for  t.  to  vary  between  t  and  t^.  Therefore, 
the  hop  or  pitch  motion  will  be  minimized  for  the  entice  counterrecoil 
stroke.  At  the  final  time  t^,  which  is  the  end  of  cointerrecoil,  the 
recoiling  parts  must  come  back  to  its  original  position  and  the  velocity 
of  the  recoiling  parts  will  be  some  specified  value  V^.  This  is  to  insure 
that  the  recoiling  parts  come  back  to  the  latch  position.  It  will  also  be 
demanded  that  the  total  cycle  time  be  equal  to  c^  seconds. 
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Formu.lating  the  above  paragraph  into  mathematical  notation  yields 


Minimize  J  ■  x  (t  } 

4  2 

subject  to  the  equality  constraints 

tj;  -  x  (t  )  -  n  +  n  =  0 
!  21  0  max 


V 2  *  x2(tf)  -  nQ  »  0 


*3  "  W  '  V£ 

ft1  «  Xjftj)  *  0 


n2  -  xjtp  «  o 


(i) 


(2) 


ft“  »  tc  -  c|T1  «  0 

i.  L 

with  the  full  set  of  initial  conditions 

x.  (0)  *  x  (0)  -  x  (0)  ■  0,  x  (0)  -  n  (3) 

*  O  4  U 

where  i  *  1,2,3  are  intermediate  and  terminal  constraint  functions  to 

be  satisfied;  ft*,  ft2,  and  ft^  define  the  times  at  which  the  intermediate 

and  terminal  constraint  functions  occur;  x^  and  x^  are  the  velocities  of 

the  recoiling  parts  and  pitch  motion  respectively;  x.  and  x  are  the  dis- 

2  4 

placements  of  the  recoiling  parts  and  pitch  motion  respectively;  i//  =  0 

is  the  constraint  on  the  displacement  of  the  recoiling  parts  such  that  at 
the  end  of  the  recoil  stroke  the  displacement  will  be  exactly  equal  to 
0  inches.  *  ■  0  is  the  constraint  demanding  that  the  recoiling  parts 

return  to  the  latch  position  at  the  end  of  counterrecoil.  •  i|>  -  0  is  the 

constraint  which  requires  that  the  velocity  of  the  recoiling  parts  come 
into  the  latch  position  at  a  velocity  in. /sec.  ft"  °  0  defines  the  time 
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at  which  the  end  of  the  recoil  occurs;  Q2  -  0  defines  the  times  at  which 
the  pitch  velocity  is  zero  and  the  one  with  the  largest  displacement  is 
selected,  thus  defining  the  time  at  which  maximum  hop  occurs;  »  0 
defines  the  total  cycle  time  to  be  exactly  equal  to  cT  seconds. 

It  was  previously  mentioned  that  in  order  to  eliminate  state  variable 
constraints  the  rod  force  was  taken  as  the  design  (control)  variable  in¬ 
stead  of  the  orifice  areas.  Using  the  rod  force  as  the  design  variable 
simplifies  the  problem  and  it  also  gives  the  engineer  more  insight  in  the 
design  process  since  he  has  an  intuitive  feel  for  the  force  levels  the 
weapon  system  he  is  designing  can  tolerate.  Thus,  immediately  the  engineer 
can  specify  an  admissible  upper  limit  for  the  rod  force  say  Rmax*  for 
his  design,  and  this  value  may  be  varied  by  the  engineer  for  any  redesign. 
The  following  inequality  constraint  must  hold  fo^  all  time  t. 

ip  *  R(t)  -  R  <  0  0  <  t  <•  t  (4) 

max  -  —  -  r 

Since  the  mathematical  model  must  represent  a  physical  re  .lization, 
to  specify  one  value  for  Rmax  is  not  enough.  This  result  was  made  avail¬ 
able  from  the  first  set  of  computer  runs  and  can  be  seen  in  Figure  4. 
Because  it  was  not  known  how  the  optimal  shape  rod  force  would  behave, 
the  design  variable  R(t)  was  allowed  to  take  on  any  shape  just  a,s  long  as 
it  did  not  exceed  R  .  Ic  can  be  seen  from  Figure  4  that  the  rod  force 
attained  its  maximum  value  at  time  tQ.  The  mathematical  model  says  that 
the  best  way  to  reduce  the  "hop"  is  to  let  the  recoiling  parts  move  for¬ 
ward  first  as  in  the  f iring-ouc-of-battery  concept.  This,  of  course,  is 
a  physical  impossibility  for  the  weapon  under  study  since  the  recoiling 
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parts  cannot  travel  forward  beyond  the  latch  position.  Additional  con¬ 
straints  were  subsequently  put  on  the  design  variable  during  the  first 
few  milliseconds  of  the  recoil  stroke. 


Figure  4.  Rod  Force  With  No  Rife.  Constraint 
'  For  First  Few  Milliseconds 

The  optimization  problem  has  now  been  formulated.  The  objective 
function  (see  Equation  II1-1)  has  been  defined  for  the  process  (see 
Equations  IV-1,2)  that  is  to  be  optimized  subject  to  the  constraints 
(see  Equations  111-2,4)  that  are  to  be  satisfied. 

All  that  must  be  done  now  is  to  put  the  problem  into  the  steepest- 
descent  formulation.  The  next  section  simplifies  the  equations  of  motion 
for  the  XM164  howitzer. 
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IV.  Translation  and  Rotational  Equations  of  Motion  for  the  XM164  Howitzer 
The  differential  equations  to  be  solved  are  given  below^Re^' 

Equations  (1)  and  (2)  are  the  translational  and  rotational  equations  of 
motion  for  the  XM164  howitzer.  Equations  (3)  and  (4)  determine  the  guide 
friction. 

M  [n  -  (4~Y  siny  +  Z  cos-,)^  “  (n  +  Y  C03y  +  Z  siny)$2] 

3  C  u  C  C 

(l) 

-  R(t)  -  B (t)  -  g  sin(yf$>)  -  u(|Sj  +  |S2|)  sgn(ii) 

{Ma(n  +  Y  cosy  t  Z^_  sin*)2  +  M^Kr^  cosy  +  -  4b  siny)2 

+  (nfe  siny  +  zt  +  4b  cosy)2]  +  Md(Y2  +  Z2)  +  1&  +  Ib  +  IdH 

+  2M  n$(n  t  Y  cosy  +  sin.)  -  M  <J>2 (n  +  Y^  cosy  +  Z  siny). 
a  t  t  at  t 

(4  -  Yfc  siny  +  Zt  cos^)  »  B(t)  ♦  (i^  -  O  +  (2) 

[R(t)  -  uOSj  +  |  S2  | )  sgn(n)  ]  (4  -  Yfc  siny  +  Z t  cosy) 

-  g  {M  [n  +  Y  cos y  t  Z.  siny]  cos(y+<p)  +  M, (Y  cosij)  -•  Z,  sin<?) 
at  t  a  a  a 

+  [Yfc  c.os4>  -  Zfc  sm<t>  +  nb  cosd+$>)  -  4b  sir. C y  +d>) ) }  -  -  c? 

[ 2r»<^  +  (n  +  Yt  cosy  +  Z^  sinyH  -  (4  +  Zfc  cosy  -  Yfc  siny)<j>2] 

(3) 

■  s  +  S  -  M_  g  COS(y-t-(J>) 

I  C  <X 


Ifl*  -  S  j  (q ,  -  r.)  +  S2  (q2  -  n)  -  B(t)  •  (4  -  4  ) 

+  R(t)  •  (4  -  42)  -  u[jSij(4  -  a)  +  |S2|(4  -  3)]  sgn(n) 


(4) 
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For  small  <f>  the  following  approximations  are  made. 
sin<j>  a  <£ 
cos4>  a  1  -  2 

The  cos(y+$)  and  the  sln(7+<p)  then  become 
cos(y+$>)  =  cosy  -  <P2  C'°f^'  ”  ^  s:*-n<' 
sin(y-H}>)  =  sinY  -  <P*  siny/2  -r  <t>  cosy. 


CONIO  ■  -  M  -g-C0N8 

ct 

C0N11  3  M„*g*C0N8/2 

d 

C0Ni2  -  M  *CON 4 
a 

CONI 3  3  -  2M 

a 

CONI  ii  *  -  ?M  •CON4 
a 


CONI 5  »  -  k*$ 

s  t 

CON16  -  H  * CON6 
a 

CON 17  -  M  *CON6*CON4 
a 

CON 18  -  -  g*M  -CON9 
a 

CON19  -  g*M  -CON9/2 
a 

CON20  -  -  g*M  *CON4*CON9 
a 

CON 21  =*  g*M  •  CON4  •  CON9 / 2 
a 

CON22  -  -  g-M.-Y, 
a  a 

CON 2 3  -  g-M  «Y,/2 
a  d 

CON 2 4  -  g*Md’Zd 

CON25  -  -  g,Mb*Yt 

CON 2 8  -  g’^-Y^ 

CON 2  7  -  g^Zf. 

CON 2  8  -  -  g’M^n^COW 

CON 2 9  -  g*Mb*nb*CON9/2 

CON30  -  g'^*;  ’CON  8 

CON31  -  -  g*Mb*?b*CON8/2 

CON 3 2  -  g*M  -CON8 
a 

CON 3 3  -  g*M  -CON4-CON8 
a 

CON34  ■  Mh'gTih'CONS 
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CON35  -  Mb*g'?o*C0N9 

CON 3 6  =>  CON 20  +  CON22  +  CON25  +  CON28  +  CON30 

CON 3 7  =  CON 21  +  CON23  +  CON26  -r  CON29  +  CON31 

CON38  »  CON 2 4  +  CON27  +  CON33  +  CON34  +  CON35  -  k 

CON 39  -  CON 15  +  CON 3 6 

C0N40  -  -  M  *g*CON9 
a 

With  the  above  definitions  Equations  (1)  and  (2)  may  now  be  written 

CONl-n  +  CON2-$  =  R(t)  -  B(t)  +  CON3  +  CONIO  t  C0N11-(J>2 
K  ^  +  Ma-n-<?2  +  CON12  *42  +  4> * CON40 - 


{M  (n  +  CON4)2  +  CONS }<|>  «  CON13-n*rr>  +  C0N14-n<f> 

+  CON38-<j>  +  CON32-n4  +  CON39  -  t-  CON16»92r, 
t  CON17*|2  +  B(t)*CON7  t  [R(t)  +  CON3)«CON6  +  CON18-n 
f  CON19*n*<{>2  f  CON37*4>2 


Equations  (5)  and  *.6)  can  be  put  ;.nto  the  following  form 


v,  ,  n  +  v  i>  ■  v 
11  12*  13 


V  4  ■  v 
229  23 


where 


v  -  CONI 


v. 


•? 


(6) 


(7) 


CON  2 
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v  «  R(t)  -  B(t)  +  CON 3  +  CONIO  +  CONll‘4.2  +  Mfln^2 

+  CON1242  +  C0N40*$ 

v  -  0 
21 

v„  «  M  (n  +  CON4)2  +  CON 5 
22  a 

V23  -  CON13*n$n  +  CON14*n$  +  CON38 •  4>  +  CON32*n<{> 

+  CON39  -  c$  +  CON16*$2n  +  C0N1?42  +  B(t)*CON7 
+  [R(t)  +  CON3]  *CON6  +  CON18*n  +  CON19*rr4>2  +  CON37-<)>2 
Equations  (7)'  can  be  written  as 
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V.  Steepest-Descent  Formulation 

The  optimal  design  problem  can  be  stated  as  follows:  Determine  the 
design  (control)  variable  R(t)  in  the  interval  0  <  t  <  t^  so  as  to 


minimize  J  “ 

MV 

subject  to 

the  constraints 

ip 

-  x  (t.) 

-  n  +  n  =  0 

i 

i  i 

0  max 

4 1 

“  X  (t£) 

1 

J 

U 

o 

2 

2  f 

0 

4>3  -  x  (tf)  -  Vf  °  0 


(2) 

1  i 

Q2  =  (t  *  0 
^  c 

*>  tc  -  c  -  o 

f  t 

$  »  R(t>  -  R  <■  0  (3) 

max  - 

and  satisfying 

x  =  t  (Equations  1V-10)  (4) 

witn  initial  conditions  » 

x.  (0)  *  x  (0)  *  x  (0)  -  0,  x  (0)  ■  r,  . 

*  3  2  0 

A.  Determination  of  the  Adjoint  Equations 

The  minimization  problem  stated  here  starts  with  an  estimated  design 
for  R(t),  analyzes  it,  and  then  improves  on  the  design.  This  steepest- 
descent  method  is  an  iterative  process  and  at  each  iteration  an  improve¬ 
ment  is  made  until  no  significant  gains  can  be  achieved.  For  a  complete 

.  .  .  . .  [Ref.  2  and  31  _  . 

development  oi  wnat  is  to  ionow,  see  '.  umy  cne  results 

of  those  derivations  will  be  used  here. 
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The  adjoint  equations  are 


0  <  t  <  t 


where  the  vectors  f  and  x  are  defined  in  Equations  IV-10  and 


X  =  t A , ,  X..,  x3,  X4]  where  T  denotes  transpose.  It  can  be  seen 

ii 


from 


Equations  (3)  that  *»  0. 
3f . 


Ill.f  [v 

3*i  V  11 


<JV 


22 


3v 


13 


3v 


v2a1[vi3'5^-  +  v22S^‘v 


23 


3v 


12 


-  v„ 


12  Sx^^  2  3  3x± 


3v 


Svn  1 


-  [v  v  -  v  V  ]  [v  -r - +  v  v  ■■  1  ]  \  /v2  v2  i  *»  1,2, 3,4 

13  22  12  23  ii  3*i  22  3xi  J  11  22 


3f 

■r — —  *  -  v  x.[C0N13*x  +  C0N14]/v 
3xt  12  l  2  11  22 

3^  “  {vnV22[2MaV13(x2  +  C0N4)  +  MaV22X3  "  v12(CON13-xiX3 


+  C0N32*x  +  C0Nl6*x2  i-  C0N18  i-  C0N19*x2)] 

<4  3  u 


-  [v  v  -  v  v  ]  [2M  v  (x  +  C0N4)]}/v2  v2 
» 3  22  12  23  a  l i  2  1 1  2 


22 


3f 

- — L  *  {v„,(2M  x  x  +  2*CON12x,)  -  v  x.  +  C0N14*x 

dx ^  22  a  t  j  3  i'  12  1 


-  C  +  2,CONl6*x  x  +  2-CON17-X  )}/v  v 

32  31122 


3f 

{v  (2*C0Nll*x  +  CON40)  -  v  (CON38  +  CON32-X 
3x2  2  4  i  ^  2 


+  2*CON19*x  x  +  2*CON37*x  )}/v,  v 

2  4  4  i  1  22 


i 
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3f  3f  3f  3f 2 

rr -2*  “  1,  t— -  ■  0,  -r —  =  0,  t —  =  0 
3x  3x,  3x  3x 

i  2  3 


3v 

23 

v  t- - v 


3v 


22 


1 


22  9^ 


23  3x 


i  J 


/V22*  1  =  1>2’3‘4> 


3f 

_ 3 

3x, 


3f3 

t— — —  x  (C0N13*x_  +  CON14) / v. 

3x.t  3'  2  2  2 

3f 

t  —  =*  lv  (C0N13*x  x  +  CON32 *x  +  CON16*x2 
3x2  22  i  3  <4  3 


+  CON18  +  CON19-X2)  -  2M  v0>(x„  +  CON4)}/v2 
^  a  23  2  22 


3 


i  2 


3  2 


3f 

-  {CON13*x.x„  +  CON14-X  -  C  +  2*CON16*xjc _+  2*CON17*xJ/v 

3f 
3x 

3f 
3x 


=  {CON38  +  CON32*x.f  2-CON19*x  x  +  2*CON37*x  }/v 


2  u 


3f 


3f 


3f 


0,  t —  13  0 ,  t —  =  1,  t —  =  0. 

3x_  3x  3x 

i.  3 


22 


22 


The  adjoint  equations  now  become 

0 


3f 

3f 

_ i 

1 

_ j 

3x 

3x 

t 

i 

!1l 

0 

3f 

3x„ 

3x 

2 

2 

3f 

3f 

i 

0 

_ 3 

3x 

3x. 

3 

3 

3f 

3f  __ 

_ 

o 

. —  ^ 

3x 

3x 

(5) 


where  the  partial  derivatives  are  defined  above. 

B.  Determination  of  the  Boundary  Conditions  for  the  Adjoint  Equations 
Because  of  the  intermediate  constraint  functions,  we  must  evaluate  X 
at  t^  ,  and  t  to  allow  for  any  discontinuities  which  may  occu'  across  t 


*F“4V***  WWBw 
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and  t  .  Since  the  initial  conditions  for  the  adjoint  equations  are  given 
at  tf,  these  equations  are  integrated  backwards..  Integration  is  carried 
out  by  integrating  from  t^  to  t2+.  Using  new  initial  conditions  at  t 
integration  is  then  performed  from  t2_  to  tJ+.  And. finally,  using  new 
initial  conditions  at  t^  integration  is  then  performed  to  tQ. 


T/MB 


Figure  5 

It  is  the  object  of  this  section  to  determine  the  initial  conditions 
at  t^,  t2_,  and  t^  for  the  four  different  integrations  performed  on  the 
adjoint  equations,  that  is,  tor  ,  <l>  and  J. 

i  2  i 

To  get  the  boundary  conditions  on  the  adjoint  equations  at  t^,  we 
choose 


•t 


z 


f 


where  f  super  and  subscripts  refer  to  the  time  at  which  the  partial 
derivatives  are  evaluated. 


For  an  arbitrary  function  Q,  we  compute 
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or 


3flf  \  •  /  9£^ 

3x  /  f  Xf  \3t  'f 


IS _ 1 

!  xf +  ( 

’  3S_\ 

^3tf) 

3Xf  j 

1  aof  \ 

1 ,!>l  \ 

k  3x  j 

!f  x‘* 1 

W)t 

a^r 

9x 


From  Equations  (2)  of  this  section  it  can  be  seen  that  0  does  not 


depend  upon  the  state  explicitly  and  therefore 


-  iS_ 


.«) 


0 .  Thus 


(tf)  3xf 


and  it  follows  from  Equations  (1)  and  (2)  that 
JT 

AJ(t  }  *  [0  0  0  0] 


<p 

X  1 


(tf) 


[0000] 


A(‘  }  -[0100] 


\lt)  -  11  0  0  01 


BOUNDARY  CONDITIONS  AT  t 


2- 


We  choose 
.  T 

X-> 


hz2  y_  ao2\ 

\3x  ~  q2  3*  / 

2~ 


where  the  superscript  2  refers  to  the  time  tr 

{iA  .  |a-  +  xT.  .  /  4  A] 


\ox  I 


\  hl  dX2  j  f 


2 
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where  a  refers  to  a  vector  of  control  parameters  and  f  for  this  problem 
is  equal  to  3,  i.e.»  we  have  t  ,  t  ,  t  ,  and  t  .  In  this  problem  there 

0  i  c  i 

are  no  control  parameters,  however,  the  additional  term  is  written  for 
completeness  of  the  expression  for  .  The  derivatives  appearing  in 
the  summation  are  zero  since  does  not  depend  upon  the  times  t^  or  t ^ . 
The  rest  of  the  terms  which  are  zero  can  be  seen  immediately  by  eval¬ 
uating  the  derivatives  in  Equations  (2).  We  how  have  that 


T 

Boundary  Conditions  for  J  »  x  (t  )  at  t,  (X ,  ®  0) 

<■>  z  i-  2+ 

„T  [0  0  0  1  ] x 

XJ-  >-[0  0  0  1] - —  [0010] 

U2-;  [0  0  1  0]x2_ 

t ,  (t,  ) 

-  [0  0  0  1]  -  [0  0  1  0] 


(7) 
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Boundary  Conditions  for  p  a  x  (t  )-n  +  n  “  0  at  t,  ( A  ■ 

i  ci  g  max  *"* 


X(J  }  -  [000  0] 


Boundary  Conditions  for  if>2  ■  x.(t^)  -  a  0  at 

,  T  .  T 

T  T  ^2  •  * 


*2  *2 
A<t2_)  “  A2+ 


*2  *2 
\t2_)  -  X(t2+) 


A2+  X2-  ~  X^2+ 

[0  0  1  01  x, 


Ip  Ip 

Y  2  •  r  2  • 

^2+  X2-  _  ^A  X^ 2+ 
<fi)(t2_) 


[0  0  1  0] 


[0010] 


Boundary  Conditions  for  iji3  «  x.t(tf)  -  •  0  at  t2_ 


A  ,1  v  -  A 
(t2_)  i-r 


,  T  T 

*3  ** 

A(t2_)  “  A(t2+) 


T  T 

V.  ip 

*2+  X2_  -  (A  *  X)2+ 

[0010]  x. 


[0010] 


A  2+  x2-  '  (A  X>2t 


[0010] 


BOUNDARY  CONDITIONS  AT  t, 


We  choose 


t  _  (  azj_  _  zj_  an1 1 
A(tj_)  \  ax  ”  ax  y 

l- 

/  a z 1 \  ao  t  /  zf  aof  \  I  zf_1  anf_1  \ 

M  ,’^r  '  l?3  ~i  2. 

|az'\  a^_  ,  t*  /  zf  anf\  zf_1  aflf-1\ 

la rl  -  aT-  (A  x)1+-  Ilf  — -  P-T  IT J 


•  /  "  -  v 


‘  I  2- 


I  az*  \  •  .  I  az‘\ 

l  ax  ]  x !  —  T  \*t  y 


In  general 


si  . Ji1  |  i  +  +  (iaii  +  isi)  +|i»i) 

3-  i.,0  \  3xl  1  8ti  /,  Is*  st  /,  V>a  / 


For  j  ■  1  we  have 


•  !  f  ao1  *  an1 1  !  an1  •  \ 

i-  “  la*  x  at  )  “  U  x  j  x_ 

If  -  +  a*' 1 -  (A1;: 

,T  3Q  ,  SI  L3X1  1+J  1_ 

X-)"8*!  1+“  ||ai.  •  j 


Boundary  Conditions  for  J  «*  x^(t2)  at  t1_ 

T  T 

jT  jT  ^+x,_  ~  (*J  x)>+ 

x(t  )  “  Ai+  "  ~i±J - : - -  t1  0  0  o] 

Ul-;  1  [1  0  0  0]x,_ 

T  T 

tT  a  .i (A  k)  , 

*<« ' »?♦  -  "ifj  ii » « oi 

*,T 

Boundary  Conditions  for  ip  ■  x  (t  )-n  +  n  «  0  at  t,  (A  *  0) 

1  2  1  0  max  1-1+ 

T  [010  0]x 

v  -  [0  1  0  0]  - - —  [1  0  0  0] 

K l-  [100  0]Xj_ 

(f2}(t  ) 

=  [0100]-  yz — \ — £ —  [1  0  0  0] 

*V(t  ) 


»,T  , 

x->"!  ■«,>(,  , 


10  0] 
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n  to 


■  *+—  >— “-'V*'1 
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T 
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Boundary  Conditions  for  i^2  =  -  nQ  a  0  at  tj_ 


ip 

i  2 

[(t  iJ 


1+ 


T  •  T* 

Al+Xl-  ~  <X  x)i+ 

[100  0]x. 

T 


[1000] 


<*!+> 

Boundary  Conditions  for  p 
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,  T 

.*3  • 
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X  f  x 
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-  (a  x) 
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[1000] 
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C.  Determination  of  the  Variation  of  the  Design  Variable 
The  variation  or  change  in  the  design  variable  R(t)  which  makes  the 
greatest  reduction  in  J,  the  hop,  is  given  by  the  following  expression 
(see^6^’  ^  an<*  where  the  desired  change  in  the  constraint  function 
is  given  by  dijj.  The  tp  constraints  of  Equations  (2)  will  in  general  not 
be  satisfied  with  the  nominal  choice  of  R(t).  Since  the  idea  is  to 


6RM  .  -  AJ(t)l 


dP2  -  di/l,  dij>  1  1/2 

_  H 

T  -I 

IT,  “  I.  ,1,  1, 

JJ  ipJ  i}m|>  ij/J 


f  W~lAV(t)l"ld4» 

u  ipy 


(15) 


drive  the  <|»  constraints  to  be  identically  equal  to  zero  along  with  mini¬ 
mizing  J,  in  the  selection  of  perturbations  the  choice  o£  the  desired  d<j/ 


will  be  -  at p.  That  is 
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■*  -af  0  <  a  <  1, 

If  a  reasonably  good  estimate  is  made  for  R(t),  the  value  of  a  may 
be  set  equal  to  1.  is  a  matrix  of  weighting  functions  whose  elements 
are  functions  of  time  which  permits  6R(t)  to  be  suppressed  in  sensitive 
regions  or  amplified  in  less  sensitive  regions.  In  this  problem  R(t) 
wa?  given  equal  weight  throughout  the  entire  recoil  and  counterrecoil 
cycle  and  was  set  equal  to  the  identity  matrix.  A  few  terms  and 
definitions  will  now  be  given  in  order  tc  evaluate  the  expression  of  6R(t). 


u^^CtjXfR.b)  3  0  0  <  t  '  t^  (21) 


0 


(22) 


Wg  is  another  weighting  matrix  and  will  be  set  equal  to  the  identity 

matrix,  b  is  a  vector  of  design  parameters  and  since  this  problem  does 

not  contain  any,  a  a  0.  Taking  variations  of  the  last  equation 

yields  the  linearized  version 

T  9j>  .  Y  9 (p  -t,  ,  T  3$  ,,  n 
p  -r2-  <Sx  +  p  “  oR  +  p  -rr~  5b  *  0. 

3x  3R  3b 

Since  0,  =»  0  we  obtain  the  following 

pT  H  6r  -  ° 

which  says  that  wherever  5R  ?  0,  p  =»  0  since  =  1,  From  (21)  it  is 

o  K 

seen  that  for  <)>  <  0,  u(t)  »  0.  However,  if  4>  is  zero  over  an  interval 
an  additional  test  must  be  satisfied. 

It  must  be  verified  that  violating  a  constraint  boundary  in  such 
an  interval  would  allow  an  improvement  in  J.  Since  $  and  R(t)  are  each 

scalars  from  (22)  we  have  that 

.T  3f 
3R 


and  it  can  be  argued  that  when  q>  =  0 ,  J  will  be  minimum  if  p  is  a  non¬ 
negative  function.  Thus,  Equations  (21)  and  (22)  provide  the  equations 

which  determine  R(t)  and  p(t).  One  more  vector,  must  be  evaluated 

oR 

now  before  6R(t)  is  determined. 
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22  3R 
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“  (v„„~  v,  *C0N6]/v  v 
9R  22  12  il  22 

3f  2  3^ 

3R~  “  °*  W  m  ° 

ill.r  !!u.  3v22  , ,  , 

3R  v22  3R  v23  3R  ^V22 

W  -  CON6/''22 

J  *± 

A  and  A  where  i  =  1,2,3  may  now  be  evaluated  by  replacing  Q  with 


J»  ipji  V  ,  and  tjy 


JL_vii;CM6\ 

V  V  1  L  V.,  3 

ii  ii  22  /  22  i 
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aQ  ^  CONfc  aQ 


1  V 
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<t>  ■  0,  p  >  0  (24) 

<J>  ■  0,  p  <  0 


I,  .,  I , ,  and  I T .  now  become 
1J1J  ijmji  JJ 

.  t ,  T  , 

I .  T  =  |  A*1'  AJdt 

-‘o 


|tf  *»TA*dt 


(25) 


ftf  JT  J 
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where  I  is  a  (3x1)  column  vector,  I  is  a  (3x3)  matrix  and  I„,  is  of 

JJ 

order  (lxl) . 
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VI.  Results  and  Conclusions 

Figur.e  (5)  represents  the  optimal  rod  force  to  minimize  hop  at  75 
degrees  quadrant  elevation  with  the  following  constraints 


R(t)  ••  22O00  lbs. 


(l; 


recoil  length  =  28  ir. 

The  resulting  hop  for  the  above  case  is  1.53  inches,  i.e.,  the  tires 
leave  the  ground  1,53  inches,  for  a  115  per  cent  maximum  rated  pressure 
breech  force.  Computer  results  indicate  that  the  present  rod  design  for 
short  recoil  yields  6.26  inches  of  hop  which  agrees  with  firing  data.  To 
obtain  the  1.53  inches  of  hop  would  require  a  redesign  of  the  orifice  areas 
for  short  recoil  and  for  counterrecoil.  One  might  question  whether  the 
resulting  curve  in  Figure  (5)  is  obtainable  with  the  XM4d  recoil  mechanism; 
if  It  is  not,  a  very  simple  solution  is  to  alter  the  curve  so  that  a  nearly 
optimal  solution  results.  If  the  constraints  were  such  that 
R(t)  "  23500  lbs. 


(2) 


recoil  length  »  29  m. 
the  resulting  hop  is  0.88  inches. 

If  one  uses  the  present  counterrecoil  groove  design  and  requires  the 
constraints  lr.  (1)  tc  hold  so  that  it  is  necessary  to  redesign  the  orifice 
areas  for  short  recoil  only,  the  resulting  hop  is  3.42  inches  or  a  45  per 
cent  reduction.  For  constraint  set  (1)  a  75  per  cent  reduction  is  achieved 
and  for  constraint  set  (2)  an  86  per  cent  reduction  results. 

The  acceleration  of  the  recoiling  parts  during  the  first  portion  of 
counterrecoil  is  an  important  factor  in  reducing  the  hop.  That  is,  the 
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faster  the  recoiling  parts  accelerate  during  this  period,  the  greater  the 
reduction  in  hop.  As  one  would  expect,  an  increase  in  recoil  length  also 
reduces  hop  significantly.  An  increase  in  the  maximum  rod  force  will  also 
reduce  hop,  for  example,  if  the  rod  force  is  allowed  to  obtain  the  value 
24160  lbs.  in  constraint  set  (1),  the  hop  can  be  reduced  an  additional 
0.32  inches.  Figure  (6)  shows  a  possible  control  rod  design  for  short 
recoil.  The  orifice  areas  were  obtained  from  the  rod  force  in  Figure  (5) • 
The  resulting  force  levels  from  the  new  groove  design  is  indicated  by  the 
dotted  lines  from  ,110  sec  to  .13  sec.  The  rod  force  is  the  same  as  the 
optimal  shaped  force  curve  from  0  to  .110  sec.  The  increase  in  hop  is 
approximately  0.1  inches.  The  recoil  length  changed  a  very  small  amount. 

An  interesting  side  point  is  that  of  the  speed  of  convergence.  The 
nominal  design  variable,  R(t),  used  for  the  first  iteration  was  such  that 
at  the  end  of  counterrecoil  the  recoiling  arts  were  250  inches  away  from 
the  latch  position  and  the  required  final  velocity  of  6  inches/sec.  was 
96  inches/sec.  In  approximately  14  iterations,  convergence  vas  obtained 
which  seems  to  be  very  fast  if  one  considers  the  complexity  of  the 
equations  involved. 

A  computational  algorithm  is  given  below. 

Step  1.  Make  an  engineering  estimate  for  R(t)  and  call  it  R°ft). 

Step  2.  Integrate  the  state  Equations  (IV-10)  with  initial  conditions 
(III-3)  and  determine  t  and  t  . 

i  i 

Step  3.  Integrate  the  adjoint  Equations  (V-5)  from  t^  to  t  with 
initial  conditions  (V-6) . 
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Step  4.  Evaluate  initial  conditions  (V-7,  8,  9,  10)  for  adjoint 
equations  at  t2  and  integrate  (V-5)  from  t  to  t 

Step  5.  Evaluate  initial  conditions  (V-ll,  12,  13,  14)  for  adjoint 
equations  at  t  and  integrate  (V-5)  from  t  to  t-. 

Step  6.  Evaluate  from  (V-24)  for  y  ,  \p  ,  p  and  J. 

i  2  i 

Step  7.  Perform  the  definite  integration  of  (V-25)  for  I  ,  1^, 
and  Ijj. 

Step  8.  Choose  dip  and  dP  in  (V-15)  where  dt^  Is  the  desired  change 
m  i (i  (V-2) . 

Step  9.  Compute  (aP)2  -  dt/  I  ‘ dif» *  If  this  quantity  is  negative, 

'P<i>  ; 

I  -•  , 1  ^ 

compute  <,  =  i(.dP)Sd^  I  dt|,J  and  replace  di^  by  r,dip. 

Step  10.  Evaluate  6R(t)  from  (V-15). 

Step  11.  Compute  new  estimate  R-(t)  ■*-  R-'(t)  +  oR(t). 

X 

Step  12-  Evaluate  gradient  squared  (Ij^,  -  I^I^l^)  for  convergence. 

Ii  near  zero,  stop;  it  not,  go  to  Step  2. 

Results  from  iirmg  tests  shew  a  significant  reduction  (50%  or  more) 
in  hop  can  be  achieved  simply  by  increasing  the  tire  pressure.  Because 
tire  performance  information  is  not  presently  available,  it  was  assumed 
throughout  this  analysis  that  the  spring  rate  or  the  tires  was  constant. 
Therefore  it  is  not  known  what  results  would  be  obtained  under  a  dynamic 
tire  response  model.  Tire  manufactures  are  looking  at  how  they  can  opti¬ 
mize  tire  characteristics  tor  the  final  configuration  in  the  tire  itself. 
In  order  to  obtain  optimum  weapon  performance  for  flexible  mount  systems 


such  information  as  tire  performance  could  be  incorporated  Into  the 
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mathematical  model  and  perhaps  Cire  characteristics  could  also  be  opti¬ 
mized  in  the  environment  for  which  they  are  being  used. 

The  technique  used  in  this  report  has  the  capability  to  optimize 
many  design  parameters  simultaneously.  If  there  exist  other  sensitive 
parameters,  consideration  should  be  given  to  optimize  them  along  with 
the  design  variable  R(t). 

This  study  clearly  indicates  that  weapon  performance  can  be  improved 
by  using  methods  of  optimal  design. 
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